A structural perspective on the dynamics of kinesin motors 

Changbong HyeorQ 

School of Computational Sciences, Korea Institute for Advanced Study, Seoul 130-722, Korea 

J. N. OnuchicQ 

Center for Theoretical Biological Physics, 
Rice University, Houston TX, 77005-1827, USA 

Abstract 

Despite significant fluctuation under thermal noise, biological machines in cells perform their 
tasks with exquisite precision. Using molecular simulation of a coarse-grained model and theoretical 
arguments we envisaged how kinesin, a prototype of biological machines, generates force and regulates 
its dynamics to sustain persistent motor action. A structure based model, which can be versatile in 
adapting its structure to external stresses while maintaining its native fold, was employed to account 
for several features of kinesin dynamics along the biochemical cycle. This analysis complements 
our current understandings of kinesin dynamics and connections to experiments. We propose a 
thermodynamic cycle for kinesin that emphasizes the mechanical and regulatory role of the neck-linker 
and clarify issues related the motor directionality, and the difference between the external stalling force 
and the internal tension responsible for the head-head coordination. The comparison between the 
thermodynamic cycle of kinesin and macroscopic heat engines highlights the importance of structural 
change as the source of work production in biomolecular machines. 
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Introduction 

In recent years, technological advances have revolutionized our understanding to biological 
systems by providing unprecedented views of biomolecular dynamics at the single molecule 
level, bringing into reality what Feynman envisioned on nanotechnology 50 years ago pQ. The 
scenery inside the cell unveiled by nanodevices displays beautiful tempo-spatial organization. 
Among a host of macromolecules that constitute the cell, of particular interest are specialized 
enzymes, namely biological motors that convert chemical energy stored in substrate molecules 
into mechanical work [2]. From transport motors that move along cytoskeletal filaments (ki- 
nesin, dynein, myosin) [3j H] to the sophisticated protein production machinery (ribosome) [51 [6] , 
biological motors play vital roles in controlling cellular functions, such as gene expression, intra 
and inter-cellular trafficking, cell motility, and mitosis [7j. A malfunction in these motors can 
cause detrimental effects in the cell. In the highly dissipative nanoscale environment that im- 
mediately randomizes any ballistic motion, it is amazing to see how these evolutionally tailored 
molecules carry out biological functions with exquisite precision. Elucidating physical principles 
that bring molecular machines into action is one of the most challenging problems in modern 
biology. 

Despite notable progress made in current nanotechnology, a precise vision on biological mo- 
tors still remains elusive due to ineluctable thermal noise in nanoscopic measurements. The 
questions begged from experiments demand additional details of the structures and dynamics 
of molecules [U |9], which require greater spatial, temporal resolutions and better control than 
current techniques can provide [10]. To this end, simulations using molecular models of bio- 
logical motors [TTHIB] can supplement the current experimental findings and make insightful 
predictions amenable to future experimental investigations. 

In the study of biological motors, however, computational approach using molecular simula- 
tions has been relatively rare compared with other area focusing on folding dynamics of small 
sized proteins. This is partly because the typical size and time scale associated with biologi- 
cal motors are far greater than usual computational approaches have dealt with. Nevertheless, 
biomolecular dynamics have a hierarchical structure in terms of their time scale, characterized 
with a large time scale separation between atom (~ fs), residues (~ ps) and domain motions 
(> fis — ms). To study the dynamics associated with a particular range of time scale, faster 
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degrees of freedom can be renormalized into an effective degree of freedom. As long as the topo- 
logical feature characterizing molecular structures are unchanged, coarse-graining the atomistic 
details does not essentially alter the global motion responsible for biological function. Thus, the 
biology can employ the strategy of condensed matter physics that simplifies the phenomenon 
of interest based on its time scale and relevant degrees of freedom [TTH |2"U] can be applied to 
the biology. Together with the current technical advances in single molecule measurements on 
diverse biological motor systems, dynamics generated from coarse-grained representation of mo- 
tors allow one to grasp the gist of design principle under which each biological motor performs 
its biological function. In this perspective, recent efforts made in studying the dynamics of 
various biological motors using coarse-grained models [m - [T8] are quite promising. 

Physical conditions that govern the dynamics of objects in the nanoscopic to microscopic 
world are fundamentally different from those for the macroscopic objects (see SI text and 
Fig. SI). To this end, the conceptual framework developed from the theory of protein fold- 
ing J2TH23] and polymer physics [25], which handle the complex dynamics of chain molecules 
over multi-scale, is of great use to decipher the design principles of biological motors. Here, 
we will discuss the dynamics of kinesin by comparing our insight gained from structure-based 
coarse grained model of kinesin with the experimental findings in each stage of biochemical 
cycle. Lastly, we will construct a thermodynamic cycle of kinesin using the neck-linker, a key 
mechanical element responsible for work generation, and highlight the importance of structural 
change in producing net mechanical work per cycle and in the regulation. 

General overview: biochemical states of kinesin 

Kinesin-1 (hereafter kinesin) transports cellular organelles along the network of cytoskeletal 
filaments and play a central role in spatially organizing the cellular environment. Two identical 
motor domains, linked via coiled-coil stalk, alternately exploit the free energy generated from 
ATP molecules to produce characteristic 8-nm steps. Single molecule experiments have revealed 
that at the physiological condition, kinesins can processively travel about 100 steps in a "hand- 
over-hand" fashion at a velocity of v ~ 1 fim/s [26ti28] . 

Recently, a growing number of kinesin structures, with various ligand states, has shed 
light on more detailed mechanisms of how kinesins function. A series of crystal struc- 
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tures suggest how kinesins adapt their conformation to varying nucleotide (NT) states via 
— >[4>— >ATP— 7-ADP-Pi— 7-ADP]— > (Fig. S2). It is currently surmised that the local conforma- 
tional changes due to the NT chemistry are amplified via the switch II helix (a4) to regulate 
(i) the neck-linker (NL) conformation and (ii) the binding affinity to MTs [29J. For kinesin-1, 
the states of NL and MT affinity are the two major components that determine the logic of bio- 
chemical cycle. Depending on the NT state, the NL is in either ordered (zippered) or disordered 
(unzippered) state; the binding affinity of kinesin head domain to MTs is either strong or weak. 
The relationship between these two binary switches in a single head, summarized in the Table [T] 
(see also Figs. S2b and S2c), forms a basic logic for the kinesin dynamics. 

The four different NT states (0, ATP, ADP -Pi, ADP) at the catalytic site enable one to assign, 
in principle, f 6 different states to the kinesin dimer; however, most of the possible dimer states 
has to be excluded. First, the ADP-ADP state (the notation (NT)£-(NT)r represents the NT 
state of trailing and leading head from left to right and we assume kinesins move from left (— end) 
to right (+ end). Below we omit (• • • )i — (• • • )r in the notations), corresponding to the "weak- 
weak" MT binding state, should be excluded since at least one of the heads should hold MTs 
tightly to remain on the MT track. Second, due to the topological constraint imposed by the NL 
zipper in the leading head (Fig. S2c), the NL of the leading head cannot be in an ordered state, 
which prohibits the 0-ATP, ATP-ATP, and ADP-ATP (or 0- ADP-Pi, ATP-ADP-Pi, and ADP- 
ADP -Pi) states from the possible dimer states. 0-ADP and <p-(f) states are also unlikely to exist 
because the time spent for ATP to dissociate from the trailing head (ATP— 7-ADP-Pi— t-ADP— > 0) 
(> I sec [30]) is longer than the timescale associated with ADP dissociation from the leading 
head (< 10 ms [30J). Therefore, only four states (ATP-ADP, ATP-0, ADP-Pi-0 and ADP-0) 
are permissible out of the 16 states. The conventionally accepted kinesin cycle is schematized 
in Fig(l|[3TH33]. 

To sustain a constant material flux along the biochemical cycle, kinesins need continuous 
supply of ATP (source) and removal of ADP and Pi (sink) (see SI text for the discussion on 
nonequilibrium steady state thermodynamics). In the kinesin cycle depicted in Fig. [TJ ATP 
binding induced stepping ([ADP-0]j_i— ?-[ATP-ADP]j) is the only mechanical motion (yellow 
arrows in Fig. [TJ that can be detected with SM measurements (the inset of Fig. [I]). All other 
steps (purple arrows in Fig. [T]) are associated with internal chemistry. In contrast to the common 
notion for macroscopic heat engines where heat released from combustion is directly used to 
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expand the volume of cylinder and provide the power stroke (see Fig. SI and SI Text), the heat 
released from the ATP hydrolysis step in the catalytic site itself is not the main driving force 
of the stepping motion. A series of structures that vary with NT state and biochemical data 
indicate that the ATP hydrolysis itself alters neither the NL state nor the affinity to MT. For 
kinesins, the free energy released from the ATP hydrolysis at the catalytic site is estimated to 
be relatively small ~ 4 k^T Even though the free energy change due to ATP hydrolysis 

in the catalytic site is small, the irreversibility of the enzymatic cycle set by the ATP hydrolysis 
has an important implication for the robust action of biological motors [37J. 

ATP binding induced conformational change 

ATP binding to a catalytic site creates new contacts to the residues that constitute the 
catalytic site (AH < 0), and reduces the local fluctuations (AS < 0). As long as the stabilization 
free energy (AG = AH — TAS) remains negative, the ATP binding is favored. Previously, 
using the molecular simulation we have shown for the catalytic domain of protein kinase A that 
local compaction of ligand binding pocket upon ATP binding induces a global open-to-closed 
transition [38]. For kinesin, despite the presence of a number of high resolution cryo-EM or 
crystal structure data [29J, [39], it is still an open question how explicitly the intramolecular signal 
transduction leads to the disorder-order transition in the NL. Nevertheless, the difference in the 
contact maps before and after the formation of NL zipper (Fig. [2^l) can be used to calculate the 
folding landscape of kinesin for each case. From the perspective of energy landscape, visualized 
in terms of the centroid of tethered head whose motion is constrained by the NL (see Fig. ^p) 
[14"] . it could be argued that ATP binding reshapes the energy landscape from the disordered 
state (Fdisorder({r})) to the ordered state (F or d er ({r})), and gives rise to a conformational force 
(fconf({r}) = — V-F({r})), which can be interpreted as the power stroke. The NL zipper formed 
in the MT-bound head biases the diffusive motion of the tethered head toward the (+)-end 
direction, which has also been termed the NL docking model \ \()\ . 

One well-known objection to the NL docking model as a driving force for the kinesin step 
is that the free energy gain calculated from EPR measurement AGdock ~ —1-2 fc#T j9j HTJ 02], 
is energetically insufficient to drive the kinesin step. However, as will be discussed in the 
next section, the energetic contribution of the NL docking (or power stroke) to the stepping 



5 



motion needs not necessarily be large as long as the NL zippered state can provide enough 
anisotropic bias, so that the tethered head can reach the next MT binding site through the 
diffusive search. Interestingly, recent SM experiments and computation have suggested that 
a structural motif called cover-neck bundle (see Fig. [2^) can guide the early stage of the NL 
docking dynamics [I3J Hlj. In addition, a simple back-of-envelope calculation using available 
thermodynamic/kinetic data provides a much greater estimate of the net stability difference 

involving the step from [ADP-</>];_i to [ATP-ADP]; as AG step « (-10 12) k B T [35l ESI SH] 

for quasi-static processes, which agrees with the work required to stall the kinesin motion W = 
(6 — 7) pN x8 nm ~ 12 ksT. AG ste p is expected to consist of AG dock and other free energy 
contributions due to the interaction between the ADP containing tethered kinesin head and 
MT surface. In fact, how tightly or loosely the free energy of molecular fuel is coupled to 
the functional motion is the key question that has been asked for many different motors [46J. 
Similarly, contrasting the mechanism of power stroke (or conformational change) with Brownian 
ratchet (or diffusive motion) may be a useful way to dissect the role of energetics and flexibility 
of molecular machines in performing their functional motion. 

Stepping dynamics 

The stepping motion in search of the next MT binding site is the most dramatic and func- 
tionally important step in the kinesin cycle. Nevertheless, because the time scale associated with 
this motion (~ (10 — 100) /xsec) is transient compared with that of the entire cycle (> 10 ms), 
details of the dynamics has been elusive and under intensive debate [HI E7J H7J EE]. The simple 
looking jump from one binding site to another, in fact, contains many complicated dynamics 
at molecular level: (i) The NL restrains the dynamics of the tethered head, (ii) Depending on 
the NT state, the NL conformation is either in disordered or ordered state, which influences the 
range of space for the tethered head to explore. The disorder-order transition of the NL in the 
MT-bound head can rectify the diffusive motion of tethered head towards the next MT binding 
site, (iii) MTs provide multiple binding sites for kinesins. Given the microscopic constraints 
and dynamics described above, what is the dynamical or structural origin of substeps observed 
in some of the experiments [271 HI]? To what portion of entire stepping motion is contributed 
by the motion due to NL docking and diffusive search? 
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Probing the stepping dynamics at molecular detail is even now an extremely demanding 
task because of the limited spatial and temporal resolutions of the current instrumentation 
[9j |27J HH 08]. Therefore, a molecular model is of great help to study this particular problem 
[!%] . Kinesin models of ordered and disordered NL, created by adjusting the strength of the 
NL native contacts (Fig. [2^,), were used to calculate the potentials of mean force (PMFs) and 
to visualize how the formation of NL zipper changes the search space for the tethered head. 
Under the disordered NL, the tethered head spans the broad area around the side-way binding 
site (c in Fig. |2]d), but cannot reach the next target binding site (e in Fig. |2}d). When the NL is 
in the docked state, the search space is confined preferentially to the forward direction relative 
to the MT-bound head, guiding the tethered head to the target binding site. A theoretical 
consideration based on two competing time scales can be proposed; if the disorder-order NL 
docking transition is slower than the time scale associated with exploring the MT surface, so that 
the tethered head can fully explore the MT surface while the NL is still in the disordered state, 
then the side-way binding site (c in Fig. [2j-b) can transiently trap the swinging tethered head; 
and consequently a substep emerges in the averaged time trace. This scenario can be simulated 
using the dynamics of a quasi-particle on a time dependent free energy surface F(r, ty/ksT = 
- log [t/r p ■ e - F °^/ fcsT + (1 - t/r p ) ■ e - F ^/ ksT ] for t < r p and F(f, t) = Fx(r) for t > t p , which 
combines the PMFs at the two extreme cases; Fo(f) under disordered NL and F\{f) under 
ordered NL. Here, t p is a parameter for the duration of NL docking transition. For a quasi- 
particle representing the centroid of tethered head with a diffusion constant D = 0.2 /im 2 /s, the 
averaged time traces exhibit the signature of a substep when r p is greater than w 20 fxs (Fig. [2^ 
left), which is similar to the one observed by Yanagida and coworkers [27] (see Fig. [2J3 right). Our 
theoretical study using multiscale simulation method suggests that the ratio between the time 
scale of NL zipper dynamics and the exploration time on MT surface determines the detailed 
pattern of the stepping [T3] , 

Interestingly, the recent experiment by Yildiz et al. for wild type kinesin using a quantum 
dot labeling [19] displayed ~ 13% of 6 nm displacement, which is likely to be the signature of 
sideway steps to the neighboring protofilament. Given the longstanding debate on the existence 
of kinesin substep [HI E7J El SB] , re-investigating kinesin time traces using an instrument with 
improved temporal and spatial resolution would be of great interst. 
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Factors determining motor directionality 



Despite a remarkable similarity among the head domain structures from 50 different ki- 
nesins belonging to the same family, the functions of kinesins are diverse; kinesin-1 moves towards 
(+)-end, kinesin-14 (C-terminal motor) moves towards (— )-end, kinesin-3 (KIF1A) undergoes 
diffusive motion along MTs [SO], and kinesin-13 quickly diffuses along MTs to reach the ends to 
depolymerize tubulins [51J. Given the structures of a motor with varying NT states and motor 
track, what determines the directionality of the motor? Brownian ratchet models suggests three 
basic conditions for the unidirectional motion: (i) asymmetric potential, (ii) thermal fluctuation, 
and (iii) athermal fluctuation (flashing ratchet or fluctuating potential) [52]. 

In kinesins, these three conditions are realized through the interplay of several factors. First, 
the interaction of kinesins in ADP-ADP state with MTs results in dissociation of an ADP from 
one of the two heads, breaking the symmetry of NT state of the two heads. Throughout the 
entire cycle, NT states remain asymmetric in the two motor heads. Second, kinesin motor 
domain binds MT in preferred direction, so that the NL in ordered state points toward the 
(+)-end of MTs. Third, as clearly demonstrated in Fig. [2Jd, the ordered state of NL restrains 
the search space of swinging head toward the (+)-end direction. Finally, the irreversibility of 
ATP hydrolysis sets the arrow of time, breaking the time reversal symmetry. 

Given the thermal and athermal (ATP-cycle) fluctuations that are provided as basic driving 
force for the conformational change in biological enzymes, the directionality of a biological motor 
is determined by the structure and conformational dynamics and its specific interaction with 
cytoskeletal filaments and other proteins [53J. 

Conformational flexibility— A strategy to avoid high free energy barriers 

A stark difference in the dynamics of biological motors from that of macroscopic machines 
is that the intrinsic energy scale responsible for the biological assembly is comparable to the 
thermal energy (see also SI Text). It is well entrenched that free energy gaps between functional 
states of biomolecules are marginal [M] , thus a biomolecule can easily adapt its conformation to 
environmental changes due to chemical transformation of molecular fuels and external tension. 
In addition, the flexibility of biomolecules makes their characteristic dynamics entirely different 
from that of macroscopic system. For instance, binding between two partner molecules such 
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as enzyme and substrate protein does not necessitate an absolute geometrical complementarity 
[55] . As suggested by Koshland, it is more likely for flexible biomolecules that interactions 
between biomolecules are realized through a mutual adaptation of the structures [56] . 

The mechanisms of partial unfolding/refolding and fly-casting to avoid large binding free 
energy barrier have been theoretically postulated [57H59J and experimentally implicated (HOI ED] 
in binding processes associated with protein-protein or DNA-protein interaction, and in proteins 
with intrinsically disordered region [62]. Our calculation suggests that kinesins also adopt this 
strategy to march along the MT by efficiently identifying the next MT binding site. Although 
the transient disruption of local structure may be difficult to detect, such dynamics could be 
important if free energy barrier associated with protein-protein recognition is too high. Indeed, 
the simulation of binding process of our kinesin model on MTs showed that the partial disruption 
of the internal structure facilitates the binding dynamics of kinesin |14j . The binding events 
monitored by using the native contacts of the kinesin's MT-binding motifs (Q p ) and the interface 
contacts between the kinesin and MT {Qi n t) show transient decrease of Q p value prior to binding 
(Fig. [3]). The transient partial local unfolding and refolding help the kinesin bind the MTs by 
bypassing otherwise a high free energy barrier. 

Out-of-phase head-head coordination in the enzymatic cycle. 

While biomolecules generally fold to function, interaction with other molecules can restrain 
a part of molecular structure, deliberately suppressing its ability to fold. Such a strategy is 
found in the design principles of biological motors that consist of multiple domains. Especially, 
dimeric kinesin maintains the asymmetrical chemical state in two motor heads and accomplish 
processive steps along MTs. 

For kinesins to maintain a high processivity, the enzymatic cycle of each head is kept out-of- 
phase. A quest to understanding the molecular origin of the head-head coordination has been 
pursued over the last decade. It has been surmised based on a series of experiments that the 
tension on NL is responsible for the out-of-phase coordination between the two motor domains. 
Nevertheless, it was not straightforwardly determined whether such coordination is realized due 
to the facilitated detachment of trailing head induced by forward tension [63] or due to the 
rearward tension induced inhibition of ATP binding to the leading head [33] 161] . Remarkably, 
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treating the MT surface as a geometrical constraint on which the kinesin dynamics occurs, 
our molecular simulation using C^-based coarse-grained model of kinesin suggests that, when 
both heads of kinesin are bound to the MT binding site as in ATP-0 or ADP-Pi-0 states, the 
NL of the leading head that is stretched backward perturbs the catalytic site of the leading 
head away from its native-like configuration, whereas the catalytic site of trailing head remains 
intact |13J . The internal tension (/ int ) exerted through the NL, estimated using the extension 
of the NL calculated from the simulation (5x) and the force-extension relationship of wormlike 
chain model (f int = k B T/l p x [1/4(1 — 5x/L) 2 — 1/4 + 8x/L]), is 8 — 15 pN depending on 
the l p value assigned |13J . The structure-function relationship implies that the deformation of 
catalytic site results in the loss of ATP binding affinity to the catalytic site, thus inhibiting the 
premature of binding of ATP. Such an ATP inhibition state is maintained as long as the two 
heads remain bound to MTs. The deformed leading head catalytic site is restored only after 
the Pi is released to change the trailing kinesin head from strong to weak affinity state with 
respect to MTs. Besides the inhibition of ATP binding, the disruption of the catalytic site in 
the leading head regulates the off-rate of ADP (k^Jj p ), which was also discussed by Uemura et 
al. [64]. Note that k^ff = 75 - 100 s" 1 [30J, 300 s" 1 [US] in the leading head but k^ff = 1 
in the trailing head. The asymmetric strain induced regulation mechanism [531 154] , the 
molecular origin of the head-head coordination, and processivity unique to the kinesin-1, is the 
consequence of the interplay between several topological constraints imposed on kinesin/MT 
complex structure. By explicitly showing the perturbed configuration at the catalytic site of the 
leading head (Fig. El), our computational study using the simple structure based model supports 
the experimental proposal of the rearward strain regulation mechanism between the two motor 
heads [131 E2 M, M EE] ■ 

Diagram of internal tension versus extension curve for the neck-linker. 

In contrast to the thermodynamic conditions imposed on macroscopic engines, molecular 
motors operate under isothermal and highly dissipative environments; thus an thermodynamic 
cycle with two distinct isotherms as in the Carnot engine is not applicable for biological motors 
(Fig. SI). Instead, conformational changes of internal structure can be linked to the thermo- 
dynamic cycle. As a simple example, the thermodynamics of the optomechanical coupling in 
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polyazopeptide |67J is best described in terms of two thermodynamically conjugate variables, 
force (/) and extension (x) (Fig. [5^,). Photon absorption/emission transforms the polymeric 
property of polyazopeptide, the persistence length (l p ) and contour length (L c ), via cis <-> trans 
conformational change. The two distinct force-extension curves with different l p and L c , cor- 
responding to the distinct conformational states of polyazopeptide, are analogous to the two 
isotherms in the Carnot cycle. The area enclosed by the optomechanical cycle is the maximal 
mechanical work that the polymer engine can extract from the photon energy when the poly- 
mer is stretched and relaxed quasi-statically. Biological machines adopt a similar strategy as 
the polyazopeptide by changing molecular conformations in response to the interactions with 
molecular fuels. In analogy to polyazopeptide whose mechanical property and conformational 
change are used to generate work, it is essential to identify relevant structural elements, respon- 
sible for internal mechanics, from its complicated architecture, in order to understand the work 
generating mechanism of each biological motor. 

For kinesin, as suggested in the previous sections, the NL is one of the key structural elements 
responsible for both the stepping dynamics and head-head regulation. Along the biochemical 
cycle, the change in motor head affinity to MTs as well as the disorder-order transition of the 
NL alters the polymer property of NL. By assuming that NL, a 15 amino acid polypeptide 
segment, is well described using worm-like chain model (see, however, Ref [68J), we propose a 
diagram in Fig. [HJd that depicts the magnitude of an internal tension (fi n t) built on one of the 
NLs (NL of red motor head) versus its extension (Sx) with its direction (5x > for forward and 
5x < for backward extension) for the full cycle. The diagram of thermodynamic cycle, built 
by tracking the state of NL along the biochemical cycle of kinesin, provides a integrated view of 
how the extension/release and conformational change of NL lead to generate work and regulate 
the head-head coordination. 

In (i)— >-(ii), the NL changes both the length and direction upon docking transition. In this 
step the residue 326 binds the /910 to form a (3 sheet, leaving only two residues at the C-terminal 
of the NL. As a result, the contour length of the NL (L c ) decreases from L c « 5.7 nm to wl.O 
nm. The L c « 1.0 nm is maintained until Pi is released from the trailing (red) head at the 
state (v). Upon ADP released at (ii)— »(iii), the leading (blue) head strongly binds the MT. In 
(iii)— ?-(iv)— >(v), Pi release from the trailing head is followed by the order-disorder transition in 
the NL, which restores the contour length into L c = 5.7 nm. If (v)— >(vi) step were to occur 
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quasi-statically, force mechanics of the NL would follow the line of the force-extension curve. 
Along (vi)— >-(vii)— >(viii) — the value of L c does not change but because the binding affinity 
of the trailing (blue) head is altered, the fi nt built on the NL is expected to vary. 

The thermodynamic cycle of kinesin by focusing on the force extension property of one of 
the NLs (Fig. [5|d) provides a few interesting points. First, the f int would be greatest when the 
two motor heads are simultaneously bound to the MT. By assuming l p w 0.6 nm for the NL, 
the fi n t for the two-head bound (2HB) state (see the section discussing the nature of ATP-wait 
state in SI text) is expected to be ff^ ~ 11 pN. It is of particular note that the estimated 
value of f1n t B ~ 11 pN should not be confused with the value of force {f stall ~ 6 — 7 pN) 
externally applied through the coiled-coil stalk that can stall the kinesin motion |49j. While 
f stall interferes with the step (i)— >(u) or (v)— ►(vi) by preventing the NL zipper process, the fl^ B 
on the NL is used to regulate the head-head coordination. In addition, the diagram suggests 
that jfnt B > f stall- Second, the diagram of thermodynamic cycle succinctly indicates the origin 
of symmetry breaking, which gives rise to the motor directionality. The finite area enclosed by 
force extension curves is formed along the process from (i) to (vi) while the process involving 
(vi)— j-(vii)— ►(viii)— >(i) creates no work. Among the series of different kinesin configurations, 
the work generating conformational change of NL occurs while the motor head is mostly in the 
trailing position (8x > 0) and its configuration is sequentially altered from ATP state to ADP 
state (see the red head in Fig. [5Jd) . Third, during the process of conformational change in red 
head spanning the half-cycle, the partner (blue) head is mostly in NT free state, waiting for the 
enzymatic cycle of the red head to be completed. It is noteworthy that although the process 
(vi)— >(vii)— >(viii) — >(i) creates no work, during which the NL of the red head points the rearward 
direction (8x < 0), the f int whose value changes between f sta n and j™ 3 is used to inhibit the 
ATP binding, which recapitulates the origin of head-head coordination. 

We used the thermodynamic cycle for the NL of kinesin (Fig. |5}d), similar to the optomechan- 
ical cycle of polyazopeptide (Fig. [5^,), to highlight the NL as a key element of work generation 
and regulation in kinesin dynamics. 
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Concluding Remarks. 

In this paper, using coarse-grained molecular simulations and theoretical considerations based 
on the theories of protein folding and polymer dynamics, we presented our perspective on how 
the kinesin, under significant thermal fluctuation, orchestrates its structure, conformational 
dynamics, and interaction with NT and MTs to constitute the characteristic biochemical cycle. 
We adopted into our coarse-grained model the hypothesis of structure-function relationship 
[69] from the study of protein function and its extension to protein-protein and protein-ligand 
interactions, and tried to address several key issues in the kinesin dynamics. Although a series 
of rate constants that constitute the enzymatic cycle may allow us to formally explain the 
behavior of kinesin motor as a function of external force or concentration of ATP [70], the 
chemical specificity of the motor itself, encoded in the microscopic rate constant and affinity to 
the MT, is determined by the structural details at each stage of the cycle. It is impressive to see 
a number of recent experimental efforts to engineer the chemical specificity of kinesin by altering 
various structural elements, which includes changing the length of NL [I9j [711 [72] , mutating the 
amino-acid at the neck region [73], and removing the negatively charged C-terminal region of 
tubulin (E-hook) [7U [7S]. The "mechanical" notions we tried to employ with our computational 
model in this paper such as stress, strain, deformation, tension, persistence length should be 
useful to clarify why the microscopic rate constant or affinity has that particular value and 
provide better idea to understand and engineer the dynamics of biological motors. 

Lastly, the notion of motor efficiency is worth mentioning. Although the literature often 
discusses the unusually high thermodynamic efficiency of biological motors [761 IH]> the ther- 
modynamic efficiency itself may not be such a relevant measure to understand the principles of 
biological motors in the cell where the amount of molecular fuels are buffered by the cellular 
metabolism. Depending on the functional goal of a biological motor, the optimization criteria 
may vary [75]. Analysis of the energy balance of kinesins [35] indicates that at physiological 
condition, 60% of the free energy stored in ATP is expended upon ATP binding to induce the 
conformational changes directly linked to the stepping dynamics; the remaining 40% of free 
energy due to the sequence of processes followed by ATP binding is partitioned into several 
microscopic steps that do not produce any work [36J. As discussed throughout the paper, the 
non-work producing steps, depicted with purple arrows in Fig. [I] or the processes corresponding 
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to (vi)— >(vii)— >(viii) — >■(!) in Fig. [5p, are responsible for regulating the MT affinity of kinesin mo- 
tor domain, the configurational state of NL, and the NT state of partner head through head-head 
coordination. Similar to the concept of kinetic proofreading to increase the fidelity of biological 
processes [37], the series of non-work producing steps, which otherwise appear to be aimless 
side reactions, are essential for kinesins to have the persistent motor function. To regulate the 
molecular configurations under excessive thermal noise, biological motors reserve substantial 
amount of free energy. Understanding the mechanism of allosteric regulation at the microscopic 
level would be one of the important topics to explore in the study of biological motors. 
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TABLE I: The NT state dependent configuration of a kinesin monomer 
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FIG. 1: Biochemical cycle of a kinesin motor. The state of kinesin changes along the cycle. The 
subscript i in the notation [NT-NT] j denotes the position on MT to which the kinesin binds. In the 
diagram the internal conformational states that do not alter the kinesin position are aligned vertically 
with purple arrows, and the step is depicted with yellow arrows. The inset shows a time trace from 
laser optical tweezers measurement. 
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FIG. 2: Stepping dynamics of kinesin. a. Difference of contact maps calculated with ordered and 
disordered NL reveals the contacts involving the NL docking and cover-neck bundle. Red (solid) arrows 
indicate the NL zipper contacts on the contact map and in the kinesin structure. Yellow (shady) arrows 
show the location of cover-neck bundle, where the N-terminal part of NL is sandwiched between two 
chain segments (residues 2-8 and 289-296). b. PMF of the tethered head on the MT surface computed 
for disordered (left) and ordered NL (right) of the MT-bound head. c. The average time traces from the 
Brownian dynamics simulation of a quasi-particle with varying Tp (the panel on the left) are compared 
with the average time traces measured using optical tweezers by Yanagida and coworkers (the three 
panels on the right) |27j . In Ref. [27] the individual SM time traces were partitioned into three groups 
depending on the stepping time ((i) t < 50 /xs , (ii) 50 fis < t < 100 /jls , (iii) t > 100 fxs) and averaged 
in each group. The panels on the right are adapted from the reference |27j . 
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FIG. 3: Facilitation of kinesin binding to MTs through partial unfolding and refolding of MT-binding 
motifs. Binding process is monitored using the fraction of native contacts within the MT-binding motifs 
made of a4, a5, a6, L12 and /35 (Q p ) and the fraction of interfacial contacts between kinesin and MT 
binding site (Qmt)- The inset shows an exemplary trajectory exhibiting partial local unfolding and 
refolding of MT-binding motifs prior to the binding. 
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FIG. 4: Internal tension regulated inhibition of ATP binding to the leading kinesin head. a. For 
the two head-bound kinesin, the internal tension (fi n t) built on the NL of the leading head disrupts 
the catalytic site and inhibits the premature binding of ATP, whereas the trailing head configuration 
is close to the native state, which is shown on the right for comparison. C^-RMS deviation of head 
domain excluding a6 helix (residues 2-315) is 1.8 A for trailing head and 3.8 A for leading head. b. 
fint value can be estimated by using force-extension relation of worm-like chain model. The extension 
of the NL in the leading head (5x) is ~ 3.1 ± 0.8 nm (distribution on the left panel). When L c = 15 
aax0.38 nm/aa~ 5.7 nm and l p = 0.5 — 1.0 nm with Sx = 3.1 nm are used, one can estimate the 
internal tension, fi nt ~ 7.5 — 12.5 pN (right panel). 
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FIG. 5: a. A molecular motor devised by using optomechanical cycle of polyazopeptide that undergoes 
cis o trans transition upon photon absorption/emission. The figure is adapted from [67]. b. The 
mechanochemical cycle of kinesins and the force versus extension curves (FECs) for the NL of red 
motor domain. There are three FECs; two for the unzippered NL (cyan) and one for the zippered NL 
(purple) . Each stage of kinesin cycle is marked with the circles (yellow) on the FECs indexed from (i) 
to (viii). 
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SUPPORTING INFORMATION 



Operational conditions for biological machines 

In the macroworld, the length, time, and energy scales of an object are well-separated from its 
surrounding, thus conventional Newtonian mechanics can determine the trajectory of the object 
deterministically. Such separation of scales, however, does not hold as the object size is reduced 
to nanoscales where the effect of thermal noise from the surrounding is substantial. While the 
thermal noise imparts incessant fluctuation to the molecule and hinders precise measurement, 
the thermal noise is an essential component to trigger the conformational changes of biological 
motors. Below we illuminate the general design principles of molecular machines by contrasting 
the operational condition of nanoscopic motors with that of macroscopic counterparts. 

Energy scales 

For biomolecules whose sizes are on the nanometer scale, the effect of thermal noise {ksT ~ 
4 x 10~ 21 J = 4pN ■ nm) on the molecule is no longer negligible because the energy scale of 
individual non-covalent bond interactions (~ 0(1) ksT) is of the same order. Dominance 
of entropy is ubiquitous in soft matters such as polymer, colloids, bubbles and biomolecules. 
Due to free energy barriers that are created by incomplete cancellation between enthalpy and 
entropy [21J, conformational transitions of biopolymers occur cooperatively. If the kinetic barrier 
associated with a transition is too large to overcome using thermal fluctuations, free energy 
borrowed from the chemical potential of molecular fuel can make the conformational change 
biologically accessible. Actively adopting the environmental noise to induce functional motion 
is one of the key design principles unique to biological motors. In contrast, for macroscopic 
motors that operate through tight couplings among a multitude of rigid bodies, the thermal 
fluctuations are rather a nuisance to achieve better accuracy and precision (IHj. An integral 
part of this review lies in understanding the soft mechanics of biological motors. 
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Thermodynamics 



Carnot engines that exemplifies the thermodynamics in the macroworld extract mechanical 
work (W) by accepting heat Qh from a hot reservoir at temperature Th and discarding the 



remaining heat Q c into a cold reservoir at T c (see Fig. SI i for gasoline engine as an example of 
practical heat engine cycle). The maximal work, which amounts to the heat transfer between the 
two reservoirs W max = Qh~ Qc, is extracted when an engine is operated quasi- statically 'without 
dissipating heat. Real heat engines, however, do not adopt the quasi-static operation for all 
practical purposes; some amount of heat dissipation is inevitable. As a result, the extracted 
work is always smaller than the heat transfer (W < Qh~Qc), thus the thermodynamic efficiency 
for the engine is bounded by r} Carnot = (Q h - Q c )/Q h = 1 - T c /T h , where Q h /Q c = T h /T c . 
To increase the efficiency, macroscopic machines are designed to minimize heat dissipation by 
lubricating the parts of the machine that are prone to produce friction. 

In contrast to the thermodynamic conditions imposed on macroscopic engines, molecular 
motors operate under isothermal and highly dissipative environments; thus an thermodynamic 
cycle with two distinct isotherms as in the Carnot engine is not applicable for biological motors. 
Instead, conformational changes of internal structure can be linked to the motor cycle (see the 
main text for the discussion using polyazopeptide and kinesin). 



Low Reynolds number environment 

Reynolds number, defined as Re = pva/rj, is equivalent to the ratio of inertial force 
{Finertiai pa 2 v 2 ) to viscous force (F viscous oc r/av), where p and 77 are the density and 
viscosity of media, and a and v are the size and speed of an object in motion |79j. Significant 
variation of Re values implies that there is fundamental difference in the dynamics. The 
Reynolds numbers associated with motions at both atomic and macroscopic scales satisfy 
Re ^> 1. A macroscopic object (a ~ cm) moving v ~ cm/s through water with the kinematic 
viscosity rj/p w 10~ 2 cm 2 j sec will have Re ~ 10 2 ^> 1; At atomic scales, the surrounding 
media looks discrete, so the dynamics such as vibrational motion of hydrogen bonds occurs 
essentially in vacuum (rj/p ~ 0), thus Re 3> 1 and inertial forces are dominant. In con- 
trast, for a typical protein, whose size is a ~ nm, the time scales associated with domain 
motions, such as looping dynamics and beta sheet formation is < (0.1 — 1) psec [80]; thus 

28 



the typical velocity of motion is v « 0.5 cm/ sec, resulting in Re ~ (10~ 5 — 1CT 6 ) <C 1. 
The estimated value of Re indicates that even for a large biomolecule like the ribosome 
(a ~ 20 nm) the molecular motions occur in a very low Re regime. In the Re <C 1 regime, 
reciprocal motions employed by a swimmer in macroworld do not gain propulsion from the 
fluid surrounding the molecule. Molecular motors, therefore, have to adopt an entirely different 
strategy from the one in the macroworld, to break time reversal and translational symmetry [79J . 

The above mentioned basic constraints enable us to draw a few general conclusions about 
molecular motors: (i) Due to constant exposure to noisy environment and structural flexibility, 
the deterministic description using the rigid body motion of macroscopic objects is no longer 
valid to describe the probabilistic nature of dynamics of biological systems, (ii) Unlike macro- 
scopic heat engines that extract work out of heat reservoirs, molecular motors transduce chemical 
energy into mechanical work by exploiting the changes of molecular topology |67J. Small and 
large structural adaptations are the vital components of free energy transduction, (iii) For 
biological motors bombarded by thermal noise, the distinction between the system of interest 
and surroundings is not obvious since system-bath interactions are non-negligible. Furthermore, 
unlike macroscopic heat engines, the heat dissipated due to friction can also be replenished by 
the thermal energy from the solvent, which is a qualitative description of fluctuation-dissipation 
theorem. Instead of minimizing the heat dissipation, molecular machines actively utilize molec- 
ular fluctuations due to thermal noise, (iv) To produce uni-directional transport or rotational 
dynamics under low Reynolds number conditions, isotropy in the diffusive dynamics has to be 
broken. 

Besides the general principles of molecular motors, understanding individual motors requires 
more specific knowledge on their architecture and related dynamics [8] . To this end, we explore 
the details of design principles in biological motors by focusing on one of the best studied 
molecular motors, kinesin-1. 
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Nonequilibrium steady state thermodynamics 

In discussing the action of molecular motors, it should be remembered that each state of the 
molecular motor along a biochemical cycle is in nonequilibrium steady state (NESS); in principle, 
one should not use the notion of equilibrium thermodynamics. Nevertheless, an extended form 
of the thermodynamics for NESS has long been suggested and recently revisited [8Tj [82]. In 
the NESS formalism, the second law of steady state thermodynamics reads TAS > Q ex or 
(W ex ) > AF where AS is the difference of Shannon entropy between the two states and Q ex 
(excess heat) is the total heat (Qtot) subtracted by the housekeeping heat (Qhk) generated in 
an infinitely slow process, where Qhk is the heat required to maintain constant non-zero flux 
between the two steady states [83J. According to the NESS thermodynamics of open systems, 
the concentrations of four NT states, held constant as a source and sinks of the system, provide 
a constant chemical potential difference that drives the biochemical cycle; one can define an 
effective free energy difference between two steady states by using the ratio of the forward and 
backward rate constants measured at near-equilibrium condition [36 | IHT1 I8~2| IM]. 

ATP- Wait state 

At high [ATP] (fa 1 mM) when the leading head is bombarded by ATP molecules every 0.5 
ms (bi-molecular rate constant for ATP binding is k° ATP = 2.0 ± 0.8/iM _1 s _1 ), the rate limiting 
step for the kinesin cycle is the hydrolysis or Pi-release from the trailing head (~ 10 ms). In 
contrast, at [ATP]w 1/iM, the leading head is inaccessible to ATP for 0.5 sec (^> 10 ms) during 
which all the processes including ATP hydrolysis, Pi release take place, thus the ATP binding 
becomes the rate limiting step. Therefore, among the series of kinesin configurations, the ADP- 
(f) state becomes the rate limiting configuration at low [ATP]. It is worth reviewing the series 
of recent studies, led by several experimental groups, to address an intriguing question as to 
whether kinesins in ADP-0 state wait for ATP binding in one head (1HB) or two head bound 
(2HB) form [HSHBH]- 

By analyzing the concentration of a/3-tubulin and kinesin-1 at various nucleotide conditions 
Alonso et al. showed that stoichiometry of kinesin-1 to tubulin is 1:1 under ADP only condition 
but this changes to 1:2 in the presence of AMP-PNP (ATP-analog) [85] . The 1:2 stoichiometry 
is consistent with the Table [I] since both NT-free and AMP-PNP states strongly bind the MT. 
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The 1:1 stoichiometry under ADP only condition results from a kinesin conformation in which 
one head is in NT-free state that binds MT strongly while the other head is still in ADP state. 
Based on the cryo-EM map, Alonso et al. proposed a structure of ATP-wait state in which the 
tethered head parks on top of its partner head, blocking its MT-binding motifs, and further 
inferred that the ATP binding to the MT-bound NT-free kinesin head unblocks the MT binding 
motifs of the tethered head, enabling the kinesin-tubulin interaction. 

Alonso et aZ.'s proposal for ATP-wait conformation (85] was examined by the more quanti- 
tative analysis using sm-FRET, which directly probed the "distance" between the two kinesin 
heads parked on MTs. In reference to the sm-FRET distributions for 2HB and 1HB struc- 
tures, the latter of which were prepared using mutant heterodimeric kinesin, Mori et al. con- 
cluded that kinesin waits ATP in 2 HB state at high [ATP] (fa ImM) and in 1 HB state at low 
[ATP](~ 2fiM). At low [ATP], FRET time trace indicates that a kinesin step has a short-lived 
2HB state, which then undergoes a transition to a long-lived 1HB state. Also, at low [ATP] the 
detached head is on average located rearward relative to the NT free MT-bound head [86] . 

The picture of ATP-wait state by Mori et al. [86] was further elaborated by Asenjo et al. 
who probed the mobility of tethered kinesin head by using fluorescence polarization microscopy 
(FPM) [87] . The FPM showed that the unbound head is highly mobile. Especially, the high 
mobility of detached head at low [ATP] rules out the parked configuration of the tethered head 
proposed by Alonso et al. In fact, the PMF calculated for the head tethered to a disordered 
NL (Fig. 2b left) succinctly demonstrates the space explored by the tethered head in 1HB state 
[14"] . The recognition of the next binding site is promoted only after the NL docking. The 
NL zippered state of the MT-bound head restricts the search space of the tethered head and 
facilitate the tethered head to locate the next MT-binding site. 

Most recently, Guydosh et al. conducted an optical tweezers experiment to provide a direct 
probe of the ATP-wait stage at low ATP concentration (=2/iM) [88] . Instead to the stalk 
region, they tethered a micro-bead to one of the two kinesin heads and probed its dynamics. In 
response to the external load that rapidly alternate between hindering (—1.7 pN) and assisting 
(+1.7 pN) phases with 14 ms time interval, the kinesin head showed signals swinging back and 
forth, which indicate that the kinesin is in 1HB state. Such a swinging motion would have 
not been possible if the kinesin had been in 2HB state. Under the assisting load (+1.7 pN) 
the tethered head showed overshoot step that reaches ~ 23 nm and recovery step (~ —7 nm) 
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upon ATP binding to the MT-bound head. The decrease of overshoot dwell with increasing 
ATP concentrations suggests that the recovery step is induced by ATP binding. Based on the 
swinging motion that freely occurs with no signature of disruption even under ±1.7 pN load, 
Guydosh et al. concluded that the ATP wait stage of kinesin is 1 HB state. The possibility of 
2 HB state was essentially precluded based on the estimate that the binding affinity of weakly 
bound ADP state is even less than thermal energy fcgT (If the location of transition state is 
assumed to be 2 nm, unbinding energy is 1.7 pN x 2 nm w 3.4 pN-nm< 4.14 pN-nm). The 
conclusion against 2HB state in Guydosh et aVs study contradicts to Mori et al. and Asenjo 
et al.'s studies [861 [87] . However, the time interval of alternating load (~ 14 ms) used in the 
experiment is still longer than the average dwell time of kinesin at saturating ATP condition 
under which 2HB state is shown to be dominant in ADP-0 state in Refs. [SSI EE]- A signature 
of disruption missing under 1.7 pN load may be due to the relatively long time interval during 
which the spontaneous dissociation of ADP containing tethered head can occur from the trailing 
head position in MTs. 
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FIG. SI: Macroscopic heat engines, a. Thermodynamic cycle of Carnot engine, b. The thermodynamic 
cycle of idealized gasoline engine (Otto cycle), which consists of (i) compression of vapor (adiabatic 
compression) — > (ii) combustion of the gasoline at the chamber with constant volume (isochoric heating) 
—> (iii) power stroke (adiabatic expansion) — > (iv) isochoric cooling. At the step (ii) heat is absorbed 
to the engine but unlike Carnot cycle the temperature involving this step is not constant. Therefore, 
the absorbed heat should be integrated over the varying temperatures. By assuming an ideal gas for 
the gasoline vapor, the maximum thermodynamic efficiency of Otto cycle is obtained with the ratio of 
volume before and after the compression, and two heat capacities (c p = (dQ/dT) p , c v = (dQ/dT) v ), 
Votto = 1 — (Vf/Vi) c " . In practice, the engine needs to be supplied with fresh gasoline vapor and 
the compression cannot be made by the single engine alone. Thus, to realize the steps (i) and (iv) the 
piston of engine is operated in concert with others. 
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FIG. S2: Kinesin structure and its conformational changes with various NT states, a. Details of kinesin 
structure. The kinesin monomer is structurally divided into three regions: the head domain (or motor 
domain) (residues 2-323), the neck-linker (NL) (residues 324-338: /39, /310) , and the neck-helix (residue 
339-: a7). With respect to the central /3-sheets ((31, (38, /33, /37, /36, /34) the head domain is divided into 
two sides; one side, made of a4, a5, a6, L8, Lll, and L12, is used to bind MTs, and the other side has 
NT-binding pocket surrounded by the structural motifs consisting of the P-loop (Nl) (86-93), switch-I 
(N2) (199-204), switch-II (N3) (232-237), and N4. Structural changes of kinesin along the biochemical 
cycle of ATP binding, hydrolysis, Pi and ADP releases are shown in b and c. b. The configurations 
of structural motifs surrounding the nucleotide binding pocket in AMPPNP (ATP analog) state (PDB 
id : li6i) (enclosed in magenta box), in ADP-AIF^ (ADP-Pi analog) state (PDB id : lvfx) (blue box), 
and in ADP-V,; (ADP analog) state (green box) |50| . The orientation of a4 (yellow helix) changes 
upon Pi release, c. Schematics showing the changes of NL configuration and kinesin-MT interaction 
indicated with red and blue arrows, respectively. 
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FIG. S3: ATP-wait state at two different ATP concentrations. At high ([ATP]=lmM) and low 
([ATP]=1 fjM.) ATP concentrations, kinesins wait for ATP in 2HB and 1HB states, respectively. At 
1HB state, the cartoon illustrate a substantial fluctuation in the tethered head [87] . 
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